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Abstract 

We show that a kilometer-scale neutrino observatory, though optimized for 
TeV to PeV energy, is sensitive to the neutrinos associated with super-EeV 
sources. These include super-heavy relics, neutrinos associated with the 
Greisen cutoff, and topological defects which are remnant cosmic structures 
associated with phase transitions in grand unified gauge theories. It is a mis- 
conception that new instruments optimized to EeV energy are required to 
do this important science, although this is not their primary goal. Because 
kilometer-scale neutrino telescopes can reject atmospheric backgrounds by 
establishing the very high energy of the signal events, they have sensitivity 
over the full solid angle, including the horizon where most of the signal is 
concentrated. This is important because up-going neutrino-induced muons, 
routinely considered in previous calculations, are absorbed by the Earth. 



PACS number(s): 95.55.Vj, 96.40.Tv, 98.70.Sa, 98.80.Cq 



I. INTRODUCTION 



It has been realized for some time that topological defects are unlikely to be the origin of 
the structure in the present Universe Therefore the observation of their decay products, 
in the form of cosmic rays or high energy neutrinos, becomes the most straightforward way 
to search for these remnant structures from grand unified phase transitions 0. Such search 
represents an example of fundamental particle physics that can only be done with cosmic 
beams. We here point out that a kilometer-scale neutrino observatory [||, such as IceCube, 
has excellent discovery potential for topological defects. The instrument can identify the 
characteristic signatures in the energy and zenith angle distribution of the signal events. It 
is a common misconception that different instruments @J^, optimized to EeV signals, are 
required to do this important science, although this is not their primary motivation. Our 
conclusions for topological defects extend to other physics associated with 10^° — 10^^ eV 
energies. 
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We will illustrate our claims by demonstrating IceCube sensitivity to: 



generic topological defects with grand-unified mass scale Mx of order 10^^ — 10^^ GeV 
and a particle decay spectrum consistent with all present observational constraints 



superheavy relics, normalized to the Z-burst scenario where the observed ultra high 
energy cosmic rays (UHECR) of ~ 10^° eV energy and above are locally produced 
by the interaction of superheavy relic neutrinos with the cosmic neutrino background 
radiation 

neutrinos produced by superheavy relics which themselves decay into the UHECRs 
and 



the flux of neutrinos produced in the interactions of UHECR cosmic rays with the 



microwave background |T3|, the so called Greisen neutrinos. This flux, which origi- 
nally inspired the concept of a kilometer-scale neutrino detector, is mostly shown for 
comparison. 

The basic reasons for our more optimistic conclusions about the sensitivity of a detector 
such as IceCube are simple. Unlike first-generation neutrino telescopes, IceCube can measure 
energy and can therefore separate very high energy signals from the low energy atmospheric 
neutrino background by energy measurement []T^ (see below). The instrument can therefore 
isolate high energy events over An solid angle, and not just in the hemisphere where the 
neutrinos are identified by their penetration of the Earth. This is of primary importance 
here because neutrinos from topological defects have energies high enough so that they 



are efficiently absorbed by the Earth [^. The signal from above and near the horizon 
typically dominates the up-going neutrino fluxes by an order of magnitude. We will show 
that the zenith angle distribution of neutrinos associated with topological defects form a 
characteristic signature for their extremely high energy origin. 



II. NEUTRINO EVENTS 

We calculate the neutrino event rates by convoluting the -|- i?^ flux from the different 
sources considered in this paper, with the probability of detecting a muon produced in a 
muon-neutrino interaction in the Earth or atmosphere: 

f f dNy 

A^events = 27r Aeflf r j J -^{E^,) Pi,^ij,{E,, , Eij_{thTesh) , COS 9^emth) dE^ d COS O^enith (1) 

where T is the observation time and ^zenith the zenith angle. We assume an effective telescope 
area of Aes = 1 km^, a conservative assumption for the very high energy neutrinos considered 
here. It is important to notice that the probability [P^^^) of detecting a muon with energy 
above a certain energy threshold i?^ (threshold), produced in a muon-neutrino interaction, 
depends on the angle of incidence of the neutrinos. This is because the distance traveled 
by a muon cannot exceed the column density of matter available for neutrino interaction, a 
condition not satisfied by very high energy neutrinos produced in the atmosphere. They are 
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absorbed by the Earth and only produce neutrinos in the ice above, or in the atmosphere 
or Earth near the horizon. The event rates in which the muon arrives at the detector 
with an energy above i?^ (threshold) =1 PeV, where the atmospheric neutrino background is 
neghgible, are shown in Table I. 

Fig. 1 shows the i^^ + fluxes used in the calculations. We flrst calculate the event rates 
corresponding to the largest flux from topological defects allowed by constraints imposed 
by the measured diffuse 7-ray background in the vicinity of 100 MeV. The corresponding 
proton flux has been normalized to the observed cosmic ray spectrum at 3 x 10^° eV; see 
Fig. 2 of reference 0. Models with p<l appear to be ruled out [[7| and hence they are not 
considered in the calculation. As an example of neutrino production by superheavy relic 
particles, we consider the model of Gelmini and Kusenko |0| . In Figs. 2 and 3 we show 
the event rates as a function of neutrino energy. We assume a muon energy threshold of 
1 PeV. We also show in both plots the event rate due to the Waxman and Bahcall bound 
T^ . This bound represents the maximal flux from astrophysical, optically thin sources, in 



which neutrinos are produced in p-p or p-7 collisions. The atmospheric neutrino events are 
not shown since they are negligible above the muon energy threshold we are using. The 
area under the curves in both Figs, is equal to the number of events for each source. In 
Fig. 4 we plot the event rates in which the produced muon arrives at the detector with an 
energy greater than i?^ (threshold). In Fig. 5 we flnally present the angular distribution of 
the neutrino events for the different very high energy neutrino sources. The characteristic 
shape of the distribution reflects the opacity of the Earth to high energy neutrinos, typically 
above ~100 TeV. The limited column density of matter in the atmosphere essentially reduces 
the rate of downgoing neutrinos to interactions in the 1.5 km of ice above the detector. The 
events are therefore concentrated near the horizontal direction corresponding to zenith angles 
close to 90°. The neutrinos predicted by the model of Gelmini and Kusenko are so energetic 
that they are even absorbed in the horizontal direction as can be seen in Fig. 5. 



Model 


^^^M+^M (downgoing) 


Nu^+u^ (upgoing) 


TD, Mx = 10^^ GeV, Qo = 6.31 x lO'^^, p=l 
TD, Mx = 10^4 GeV, Qo = 6.31 x 10-^5, p=2 
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1 

0.3 


TD, Mx = 10^5 GeV, Qo = 1.58 x 10-^4^ p=l 
TD, Mx = 10^5 GeV, Qo = 1.12 x lO^^^^ p=2 


9 
2 


1 

0.2 


Superheavy Relics Gelmini et al. [|10 


30 


1.5 X 10-^ 


Superheavy Relics Berezinsky et al. |11| 


2 


0.2 


Superheavy Relics Birkel et al. |12| 


1.5 


0.3 


P-7CMB (^^max = 2.2) 13| 


1.5 


1.2 X 10-2 


W-B limit 2 x 10"^ E'"^ (cm^ s sr GeV)"^ 


8.5 


2 


Atmospheric background 


2.4 X 10^^ 


1.3 X 10^2 



Table I: Neutrino event rates (per year per km^ in 27r sr) in which the produced muon arrives 
at the detector with an energy above i?^ (threshold) =1 PeV. Different neutrino sources have 
been considered. The topological defect models (TD) correspond to highest injection rates 
Qo (ergs cm"^ s"^) allowed in Fig. 2 of 0. Also shown is the number of events from p-7cMB 
interactions in which protons are propagated up to a maximum redshift Zmax = 2.2 |T^ 



and the number of neutrinos from the Waxman and Bahcall limit on the diffuse flux from 
optically thin sources |]16[- The number of atmospheric background events above 1 PeV 
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is also shown. The second column corresponds to downward going neutrinos (in 27r sr). 
The third column gives the number of upward going events (in 27r sr). We have taken into 
account absorption in the Earth according to reference . IceCube will detect the sum of 
the event rates given in the last two columns. 

Energy measurement is critical for achieving the sensitivity of the detector claimed. For 
muons, the energy resolution of IceCube is anticipated to be 25% in the logarithm of the 
energy, possibly better. The detector is able to determine energy to better than an order 
of magnitude, sufficient for the separation of EeV signals from atmospheric neutrinos with 
energies below 100 TeV. Notice that one should also be able to identify electromagnetic 
showers initiated by electron and tau-neutrinos. Their energy measurement is linear and 
expected to be better than 20%. Such EeV events will be gold-plated, unfortunately their 
fluxes are expected to be even lower. For instance for the first TD model in Table I (p=l, 
Mx = 10^^ GeV and Qo = 6.31 x 10~^^ ergs cm~^ s~^), we expect ~ 1 contained shower 
per year per km^ above 1 PeV initiated in charged current interactions of z/g + v^- The 
corresponding number for the Gelmini and Kusenko flux is ~ 4 yr~^ km~^. 

One should also worry about the fact that a very high energy muon may enter the de- 
tector with reduced energy because of energy losses. It could become indistinguishable from 
atmospheric background [T^. We have accounted for the ionization as well as catastrophic 
muon energy losses which are incorporated in the calculation of the range of the muon. In 
the PeV regime region this energy reduction is roughly one order of magnitude, it should be 
less for the higher energies considered here. 

In conclusion, if the fluxes predicted by our sample of models for neutrino production in 
the super-EeV region are representative, they should be revealed by the IceCube observatory 
operated over several years. 
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FIG. 1. Maximal predictions of + v'fj, fluxes from topological defect models by Protheroe and 
Stanev (p=l,2). Also shown is the z^^ + i^^ from superheavy relic particles by Gelmini and Kusenko 
and the flux by Berezinsky et al. 
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FIG. 2. Differential + event rates in IceCube from the topological defect fluxes in Fig.l. 
The muon threshold is i?^ (threshold) =1 PeV. We have separated the contribution from upgoing 
and downgoing events to stress the different behavior with energy. The event rate expected from 
the Waxman and Bahcall bound (see text) is also shown for illustrative purposes. The rate due to 
atmospheric neutrinos is negligible (see Table I) and hence it is not plotted. 
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FIG. 3. Differential V/j, + event rates in IceCube from super-heavy relic particles. We have 

separated the contribution from upgoing and downgoing events to stress the different behavior 
with energy. The muon threshold is i?^ (threshold) = 1 PeV. The event rate due to atmospheric 
neutrinos as well as the one expected from the Waxman and Bahcall bound (see text) is shown for 
illustrative purposes. The rate due to atmospheric neutrinos is negligible (see Table I) and hence 
it is not plotted. 
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FIG. 4. i/^ + i>ij, event rates in IceCube from the fluxes in Fig.l. The plot shows the number 
of events in which the produced muon arrives at the detector with an energy above i?^ (thresh). 
Atmospheric neutrino events and the event rate expected from the Waxman and Bahcall upper 
bound (see text) are also plotted. The topological defect (TD) models shown (p=l and p=2) 
correspond to Mx = 10^^ GeV. Upgoing and downgoing events are shown separately. 
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FIG. 5. Zenith angle distribution of the + Ufj, events in IceCube in which the produced 
muon arrives at the detector with energy above 1 PeV. Left: Topological defect models. Right: 
Superheavy relics, cos(^zcnith) = ^1 corresponds to vertical upgoing neutrinos, cos(^zcnith) = to 
horizontal neutrinos and cos (^zenith) = 1 to vertical downgoing neutrinos. The detector is located 
at a depth of 1.8 km in the ice. 
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